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I .  Introduction 

The  understanding  of  elementary  excitation  spectrums  of  two-  and  three- 

4 

dimensional  liquid  He  has  been  one  of  the  important  problems  in  low- 

temperature  condensed  matter  physics.  After  Landau's  well-known 

phenomenological  theory^  for  liquid  ^He,  many  microscopic  theories  have  been 

developed  by  Bogoliubov,  Feynman  and  Cohen,  Lee  and  Yang  and  others.  By 

considering  of  interactions  between  quasiparticles,  Kebukawa  et  al  and 

Iwamoto  have  tried  to  explain  the  multi-excitation  spectrum  in  bulk  liquid 
4 

He.  However,  most  theories  have  been  primarily  developed  for  absolute  zero 
temperature  and  have  yielded  the  normal  dispersion  (convex  down)  relation, 
which  yields  the  wrong  spectrum,  whereas  the  anomalous  dispersion  (concave  up) 
relation  is  correct.  In  the  former  case,  four-phonon  processes  turn  out  to  be 
the  lowest  order,  while  in  the  latter  case  three -phonon  processes  are  the 
lowest. ^ 

8 

Recently  Chin  and  Krotscheck  computed  the  ground- state  structure  and 

.  ,  4 

collective  excitation  energies  of  He  droplets  at  zero  temperature,  which  are 

q 

described  by  a  generalized  Feynman  theory,  while  Krishna  and  Whaley 

calculated  the  excitation  spectra  of  compressional  modes  of  ^He  for  20,  70  and 

240  clusters  at  0  K.  Chakraboty  et  al^*^  derived  the  excitation  spectrum  from 

the  random-phase  approximation  with  the  assumption  that  Bogoliubov  excitation 

is  noninteracting,  and  they  found  that  the  spectrum  agrees  only  qualitatively 

with  the  experiment.  On  the  other  hand,  Reppy  et  al  observed  sharp  cusplike 

heat  capacity  singularities  in  coincidence  with  the  superfluid  transition, 

which  is  strongly  related  to  the  excitation  spectrum,  and  confirmed  that  ^He- 

filled  aerogels  are  not  in  the  same  universal  class  as  bulk  helium.  Stirling 
13 

et  al  performed  high-precision  neutron  scattering  measurements  of  the 
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temperature  dependence  of  the  phonon  and  roton  excitation  in  liquid  He  at 
saturated  vapor  pressure  in  both  the  superfluid  and  normal -fluid  phases. 

For  two-dimensional  liquid  ^He,  in  recent  developments  of  experimental 

4 

techniques  ,  the  properties  of  He  films  have  been  widely  investigated.  Some 
years  ago,  Isihara  and  Urn  derived  successfully  the  elementary  excitation 
spectrums,  that  are  microscopic  only  in  the  long  wavelength  limit,  for  two- 
and  three-dimensional  liquid  ^He  in  the  zero  temperature  limit  using  the  ring 
diagram  approximation,  which  gives  anomalous  phonon- like  behavior  at  low 
momenta  and  roton-like  at  high  momenta.  Starting  with  these  elementary 

excitation,  we  have  successfully  explained  phonon  decay, first, second 

17  18  1 9 

and  third  sound,  their  sound  attenuation,  thermal  conductivity,  first 

viscosity^^  and  thermal  diffusion. 

22  23 

According  to  the  experiments  performed  by  Wood  and  others,  it  is 
well  known  that  the  temperature  variation  of  the  elementary  excitation  is 


small  but  becomes  significant  in  the  vicinity  of  the  A-point.  However,  we 
have  not  found  any  calculations  of  this  temperature  variation  in  two- 
dimensional  liquid  ^He.  All  the  more,  we  still  do  not  know  whether  or  not  the 
effect  of  the  temperature  variation  on  three-phonon  processes  is  significant. 


which  is  a  question  not  addressed  by  Landau  or  Khalatnikov. 

It  is  the  purpose  of  this  paper  to  evaluate  the  temperature  variation  of 
the  elementary  excitation  spectrxan  of  two-dimensional  liquid  ^He  on  the  basis 
of  a  microscopic  theory  in  the  long  wavelength  limit  in  the  ring  diagram 


approximation,  and  then  to  investigate  the  temperature  effect  on  the  structure 
factor,  thermodynamic  functions,  fluid  density  and  various  sounds  of  two- 
dimensional  liquid  ^He.  In  Sec.  II  we  shall  treat  the  pair  distribution 
function  using  the  chain  diagram  approximation  and  present  some  basic  formulas 
which  will  be  used  in  the  following  sections.  The  structure  factor  and  the 


elementary  excitation  spectrum  will  be  evaluated  in  Secs.  Ill  and  IV, 
respectively.  The  thermodynamic  functions  and  fluid  densities  will  be  given 
in  Sec.  V.  Ue  shall  evaluate  the  first,  second  and  third  sounds  in  Sec.  VI, 
and  finally  in  Sec.  VII  we  shall  discuss  our  present  results  in  comparison 
with  related  work. 


II.  Chain  diagram  approximation 

In  this  section  we  offer  the  formula  of  the  pair  distribution  function 
within  the  chain  diagram  approximation  and  related  formulas  and  then  derive 
Che  excitation  spectrum  in  the  long  wavelength  limit  at  absolute  zero 
temperature.  The  chain  diagram  approximation  is  justified  since  in  a  two- 
dimensional  Bose  system  the  pair  distribution  function  depends  on  the  modified 
Bessel  function  of  order  zero,  i.e.,  K^Cr),  in  contrast  to  an  ideal  gas. 
Within  this  approximation,^^  the  pair  distribution  function  ® 

dimensional  Bose  liquid  is  given  as 

p-(r)  -  n^  +  I,(r)  -  - 

where  n  is  the  number  density,  P  -  1/k  T,  u(q)  is  the  Fourier  transform  of  the 

O 

interaction  potential,  is  the  j-th  eigenvalue  of  the  effective  boson 
propagator  representing  the  unit  of  a  chain,  and  l2(^)  the  ideal-gas 
contribution  given  as 


f  (q) 

j  1  +  u(q)A  (q) 


(2.1) 


(2.2) 
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The  eigenvalues  (q)  of  the  effective  propagator  can  be  obtained  from  the 
expression 


A,(q)  - 


da  f  dp  f(p)[l+f(p+q)]exp[a{p^-(p+q)^)]e  ^  ,  (2.3) 


where  f(p)  is  the  Bose-Einstein  distribution  function 


f(p)  - 


1  -  ze 


where  z  is  the  fugacity. 

The  integrand  of  Eq.  (2.1)  can  be  divided  into  two  parts  as 


(2.4) 


/  1  +  u(q)A.(q)  "  ■  2  1  +  u(q)A  (q) 


A,(q) 


(2.5) 


From  the  definition  of  Aj(q)  in  Eq.  (2.3),  the  first  term  on  the  right-hand 
side  of  Eq.  (2.5)  yields 


(2.6) 


Hence,  the  pair  distribution  function  becomes 


1  *  u<,)Aj 


P^{t)  -  n  -  n6 


(2.7) 
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The  ideal-gas  term  is  eliminated,  and  it  is  much  easier  to  treat  Eq .  (2.7) 
than  Eq.  (2.1)  because  the  numerator  of  integrand  in  the  latter  is  linear  in 

Aj(q). 

In  terms  of  the  pair  distribution  function,  the  internal  energy  is  given 
by 


U(T)  -  U  (T)  +  f 
o  2 


dr  ^(r) 


2_ 

3/9 


(2.8) 


where  ^  is  a  coupling  constant,  U^(T)  is  the  ideal-gas  energy,  A  is  the  two- 
dimensional  volume,  and  (^  ( r )  is  the  two-body  interaction  potential. 
Substituting  Eq.  (2.7)  for  the  pair  distribution  function  into  Eq .  (2.8),  we 
obtain  the  internal  energy  as 


U(T)  -  U  (T) 
o 


vlh 


I 


dr  ^(r) 


nA 


J 


-da. 


(2ir)- 


u(q) 


+ 


A 

2  3/3 


-da _ 


u(q)-^^(q)  - 

1  +  ?u(q)Aj (q) 


(2.9) 


The  effective  eigenvalues  (q) 
given  by 


evaluated  to  first-order 


from  Eq . 


(2.3)  are 


Aj(q) 


2 

q 


2  no 


2  ’ 


(2.10) 


and  substitution  of  Eq.  (2.10)  into  (2.9)  yields  the  final  expression  for  the 
internal  energy  at  absolute  zero  temperature  as 
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U(T)  -  ^  n^A  I  dr  ^(r)  +  f  f  (£(q)  -  -  nu(q)  ] 

J  J  (2n) 

+  A  f  £(q)f(E)  .  (2.11) 

^  (27t)^ 

BE 

where  f(£)  -  l/[e^  -  1]  is  the  Bose  distribution  function,  and 

^(q)  “  [q^  +  2nu(q)q^]^  (2,12) 

is  the  excitation  energy  which  was  given  originally  by  Bogoliubov  and 
Zubarev. 

In  the  expression  for  the  internal  energy  given  by  Eq.  (2.11),  the  first 
and  second  terms  represent  the  ground-state  energy  and  the  quasiparticle 
excitation  energy,  where  the  latter  demonstrates  a  Landau-type  excitaiton  but 
does  not  provide  for  the  temperature  variation  of  the  elementary  excitation 
spectrum.  For  finite  temperatures,  the  above  approach  is  not  quite  correct; 
the  eigenvalues  must  have  terms  which  are  dependent  on  the  de  Broglie  thermal 
wavelength,  and  we  shall  discuss  this  problem  in  the  following  sections. 

III.  Structure  factor 

In  the  chain  diagram  approximation,  the  structure  factor,  S(q),  is  given 
oy  the  integrand  of  Eq.  (2.7)  as 

’  ^n  )  1  +  u(q)A  (q)  ’ 

j  ^ 

Combining  Eqs.  (2.10)  and  (3.1),  we  obtain 


8 


q  coth(^)3q)  [  [q^  +  2nu(q)]^ 

S(q)  -  - - T -  .  (3.2) 

[q  +  2nu(q)] 

26 

This  result  agrees  with  those  obtained  by  others.  In  order  to  extend  this 
calculation  to  higher  temeperatures ,  we  may  use 

Aj(q)  -  A°(q)  +  AAj(q)  .  (3.3) 


where 


Xj(q)  - 


2  no 


12n-Y(z)q 


fq"- 


(3.4) 


and 


AAj (q) 


ISirCz^nQ' 


^[q 


2x1.2, 


7(z) 


62(2) 

C^(z) 


0^(2) 


(3.5) 


The  structure  factor  to  first  order  in  AA  is 


JL  \  _ 

^  1  +  u(q)A°(q)  ^  [1  +  u(q)Aj(q)]^ 


(3.6) 


Substituting  Eqs .  (3.  A) -(3. 5)  into  (3.6)  and  summing  over  j,  we  obtain  (see 


Appendix  A) 
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S(q) 


^q  coth(-;3q)[q  +  ^ 

..  2nu^a)  +0!  1 

^  2[q2  +  1 

2[(nu(q))^  +  nu(q)d]^-^ 

q  coth(^;3q)  [q^  +  U  J^  , 

2fq^  +  U 


- 2nu(q) 

2[(nu(q))^  + 


nu(q)a] 


{: 


gq^  coth(^^q)[q^  +  UJ^ 


12[(nu(q))^  +  nu(q)a]‘^[q^  +  U  ] 


,3/2 


3  1  2  ^ 

_ aq  coth(-;9q)[q  + 

12[(nu(q))^  +  nu(q)Q]‘^[q^  + 

r  aq^(^)^csech^(^/9q)[q^  4-  U 
''6/9[(nu(q))^  +  nu(q)al^[q^  +  u  ] 

4  2  1  2 

gq  csech  (j;3q)[q  +  U^]h 

6^[(nu(q))^  +■  nu(q)a]^[q^  +  u^]/ 


(3.7) 


where 


a(0,z)  - 

P 

-  nu(q)(l  +  [1  +  g(^,2)/nu(q)]^) 
U_  -  nu(q)(l  -  [1  +  a(/5,2)/nu(q)]‘') 
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In  Eq.  (3.7),  the  first  bracket  represents  the  first-order  modifications,  and 

the  second  bracket  shows  the  AA .  correction.  We  will  later  discuss  these 

J 

terms  in  detail  in  Sec.  VII. 

IV.  Elementary  excitation  spectrum 

In  Sec.  II  we  used  Eq .  (2.10),  evaluated  in  the  first-order 

approximation,  to  obtain  the  excitation  spectrum  at  absolute  zero  temperature. 
In  this  section  we  evaluate  the  temperature  variation  of  the  excitation 
spetcrum.  Within  the  chain  diagram  approximation  for  the  pair  distribution 
function,  we  have  the  following  expression  for  internal  energy: 

U(/9)  -  U  (^)  +  ^  fr  [^n^  f  dr  ,^(r)  -  y9n  f  ,  u(q)] 

o  2  J  J  (2x)^ 

The  above  equation  for  the  internal  energy  is  equivalent  to  the  result  of  the 

ring  diagram.  Making  use  of  the  first-order  approximation  for  A^  and  summing 

over  j ,  we  arrive  at 


U(^)  -  U^(^)  +  ^  J 


dr  ^(r)  +•  ^nn'A  J  dr  <^(r) 


I 


+  f  coth(^)  +  nu(q)  +  /9n’u(q) 

J  (2jr)  '■ 


-  q[q^  +  2nu(q)]^  coth  x  (q^  +  2nu(q)]^ 
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2 

^nu(q)q  coth(^^)[q^  +  2nu(q)]^ 
[q^  +  2ru(q)]^ 


}  • 


(4.2) 


where  n'  -  [dn/dfi] 

The  expressions  of  interest  in  Eq.  (4.2)  are  the  terms  that  have  the 
form  2^(q)  coth[£(q)/2] ,  which  can  be  reexpressed  as 

^E(q)  +  E(q)f(E)  ,  (4.3) 

where  £(E)  is  the  Bose  distribution  function.  As  can  be  readily  confirmed, 
the  total  energy  in  Eq.  (4.2)  consists  of  two  parts:  one  corresponding  to  the 
ground-state  energy,  i.e. ,  the  zero-point  energy,  and  the  other  to 
quasiparticle  excitation  energy.  Since  we  are  mainly  interested  in  the 
quasiparticle  excitation  energy,  we  will  discard  the  constant  ground-state 
energy  term. 

Extending  this  calculation  to  higher-order  in  the  temperature  by 
including  terms  of  order  ^  ^  in  the  eigenvalues,  and  using  Eq.  (3.3)  together 
with  Eqs.  (3.4)  and  (3.5)  and  summing  over  j,  we  obtain 

2  ^ 

U(^)  -  U  (^)  +  ^  f  d?  ^(r)  +  ^nn'Al  d^  ^(r)  -  f  |  nu(q) 

o  2  J  J  ^  J 


f  dq  „  ,  .  .  ^  A  f  dq  f _ nu(q) 

2  ^^^^  ^2!  21  2 

J  (2x)^  ^  (2n)  '■[(nu(q))  +  ni 


nu(q)a 


2  H 

4  qfq  + 


2  12 

-  ^  U  nu(q)q  a  coth(-^q  ) 

X  coth(^)[q^  +  U  -  - ; - 

[(nu(q)  +  nu(q)a] 


12 


nu(Q')a 


coth(fr [q^  +  U^] 


[(nu(q))^  +  nu(q)a]^  [q^  + 


nu(q)aq~ 


coth(f^)[q^  ^+1^ 


[(nu(q))^  +  nu(q)a)\j^  [q^  +  U^] 


/3q  coth(~)  [q 


fq^  +  U 


2[q‘-  +  U^l 


'  +'  fn'u(q)  +  _ +  n*u(q)a  +  nu(a)oO| 

+  U  ^  2f(nu(q))^  +  nu(q)a]^ 


A  r  .  r 
2  J  ,o,^2  r 


nu(q)q[q^  +  U 


(2it)^  ^  [(nu(q))^  +  nu(q)a] 


-  coth(|^)[q^  + 


2  12 
nu(q)aq  coth(j  ^q  ) 


nu(q)aq 


[(nu(q))^  +  nu(q)a]^  l(nu(q))^  +  nu(q)a]^ 


coth(^)[q^  +  L_ 
[q^  +  U  l"" 


nu(q)aq~ 


cothcfr  [q^  +  UJ^ 


[(nuvq))^  +  nu(q)a]^  [q^  +  U 


^q  coth(^)[q^  +  U  l""  . 

+  2 - 4  h'u(q) 

2[q^  +  U 


2nn'u(q)  +  u'u(a)a  +  nufa)o 
2[n^u(q)^  +  nu(q)a]^ 


^)}  ■ 


(4.4) 


where  a'  -  (da/dB)  As  we  did  before,  taking  only  the  qvias iparticle 

z 

excitation  spectrum  from  Eq.  (4.4),  we  get 


E(T)  - 


j  nu(q)q[q  +  U^] 

M (nu(q) )^  +  nu(q)a]^ 


nu(q)aq 


[(nu(q))^  +  nu(q)a]^^ 
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nu(Q)Qq 


nu(q)aq' 


[(nu(q))^  +  nu(q)a]^[q^  +  [(nu(q))^  +■  nu(q)a]^[q^  + 


_ Sn  ,  ,  Lj^,  ^  2nn*u(q)  +  n'u(q)a  +  Tiu(q)a'l| 

2[q^  +  ^  2t(nu(q))^  +  nu(q)a]^ 


nu(q)  q[q^  +  U 


r  nu(,q;  qiq 
^  [(nu(q))^  + 


nu(q)aq 


nu(q)a]^  [(nu(q))^  +  nu(q)a]^ 


,nuCa)Qq„ 


nu(q)aq' 


[(nu(q))^  +  nu(q)a]^[q^  +  U  [(nu(q))^  +  nu(q)a]^[q^  +  U 


.  _ Sq  f  ,  2nn*u(q)  +  n'u(q)a  -t-  nu(q)aO| 

2  W  2  h  I  ■ 

2[q  +  U  1^  2[{nu(q))^  +  nu(q)a]^ 


(^.5) 


In  Eq.  (4.3),  Che  terms  inside  the  first  bracket  represent  the  temperature- 
dependent  phonon  spectrum,  while  the  second  bracket  corresponds  to  the  free 
particle  in  nature.  Equation  (4.5)  reduces  to  the  excitation  spectrum  [Eq. 
(2.12)]  at  absolute  zero  temperature  because  a()3,z)  tends  to  zero  in  this 
limit.  Since  a  depends  on  Bose  statistics,  the  excitation  spectrum  depends  on 
statistics . 

To  have  a  more  explicit  form  for  the  excitation  spectrum,  we  may  choose 
a  mock  potential.  However,  we  have  adopted  a  soft  potential  with  a  Lennard- 
Jones-type  tail  in  our  previous  works,  which  we  introduce  again  here,  i.e.. 


4(r)  - 


r  <  a 


-  (J)‘]  , 


r  >  a 


(4.6) 
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From  the  Fourier  transform  of  this  soft  potential,  we  may  find  an  expression^^ 
for  small  q , 


2nu(q)  -  A^(0)  +  A^q^  +  A^q^  -  A^q^ln(qa)  -  A  q^  +  . 


(4.7) 


where  the  coefficients  A^^  are  given  as 


a2(0)  -  2xna2(V.  ^  h  EJ  ,  ~  ^  (h  -  V  )  , 

J-  z  2  o  o 


o  10  ^o' 


a2  ,  6 

A2  —  4xna 


■  10  x\o.,  X  10^)}  ; 


(4.8) 

7  -  0.5771 


For  small  q  and  the  low- temperature  region,  substituting  Eq.  (4.7)  into 

•  (4.5),  the  temperature  variation  of  the  phonon  spectrum  can  be  evaluated 
from  the  expansion 


f(q,T)  -  (A^  +  ^)  q  -  q2  +  (^2  ^  ^^^^3 


8A“ 


A^ 

o 


-  —  f 

8a'  I 


a2  a"^ 

-2  ^1. 


4A"  2A 


^]] 


5 

q  + 


(4.9) 


where 


B  -  (1  +  A,)/2A 
i  o 


.0-^2  lil^l 
2*  L  2  ‘A„  J  ■ 
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We  can  easily  confirm  that  Eq.  (4.9)  is  the  elementary  excitation  spectrum  of 

the  Bose  liquid  at  very  low  momenta,  and  as  T  -»  0 ,  a  tends  to  zero,  so  that 

Eq .  (4.9)  reduces  to  what  we  have  derived  for  the  spectrum  at  absolute  zero 
14 

temperature , 


^(q)  -  A(0)q  +  Bq^  +  Cq^  +  ...  . 


(4.10) 


For  large  q,  we  take  only  the  dominant  term  to  get 


£■(9)  -  q[q^  +  (4jrna^/qa)V^J^(qa) 


(4.11) 


The  Bessel  function  is  oscillatory,  and  hence  there  can  be  a  minimum  in  the 
energy  curve.  Around  this  minimum  one  can  find  the  familiar  roton  spectrum^ 


2 


(4.12) 


where  m*,  A  and  q^  are  the  roton  effective  mass,  energy  gap  and  minimum  point, 
respectively. 


V.  Thermodynamic  functions 

In  the  view  of  the  energy  spectrum,  we  follow  Landau's  theory  to  obtain 
the  phonon  and  roton  energies.  Using  the  small  and  large  q  expansions  of  the 
energies,  i.e.,  Eqs.  (4.9)  and  (4.12),  we  obtain  the  thermodynamic  functions 
in  a  straightforward  fashion  as 


A^O)  “ 


a‘(0) 


<v> 
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A^(0) 


[1  + 


3r(A)>r(Z) 

4r(5)B 


Kkji)^ 


+ 


^!CLZ}g 

A^O) 


<v>' 


E  n) 

rot 


^*k  T.4 

-r-^] 


1  2 

qo[A  +  ^  (kgT)"]  e 


The  corresponding  specific  heats  are  given  by 


S  <’■>  -  2^  P 

ph  ^  A‘(0) 


^  (  2  ^  2.41aIZ^3I  ^  3 

B  B 


4‘5.!C(5)P 

A^(0) 


AAuT 


SoP’  ■  •'b-oK  *  -  " 


V  V 


The  theoretical  expressions  for  other  thennodynajnical  functions  are 


2>r  L2 


a2(0)  ®  A^0)  ® 


4!.f(5)P  r  37CZK(41 
a5(0)  ^  ■  A(0)Br(5) 


1“."‘ ■  XT'*-"*  •  ■ ) 


F  ^(T) 
rot 


(T) 


2x 


rami 

V(0) 


<kjT)' 


4!7(?)C..a.} 

A^(0) 


(kgT) 
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.  (5.1) 


(5.2) 


■}  • 
(5.3) 


(5.4) 


.  (5.5) 


(5.6) 


5£im  fi  4 

L  ■  BA(0)f(5)J^V^ 


A^(0) 


■<v>' 


A^O)  ® 


■}  • 


(5.7) 


-  0^)“  ■‘B-eg 


"  Vi' 


-A/kgT 


(5.8) 


{^^(kgT) 

^a2(0)  ^ 


3  . 

A^(0)  “ 


2*4!r(5)B, 


A^(0) 


L"-  A(0)Bj^f(5)j^‘'B^^  •■•  J  ’ 


(5.9) 


■aAbT 


(5.10) 


If  the  quas iparcicles  move  with  an  average  drift  velocity  with  respect 
to  the  rest  system,  and  this  momentum  is  associated  with  the  normal  fluid, 
then  its  density  is  given  by 


WTl  .  -  iw  r  .3,  f3f(E(a)l1 

J„  ^  3E<q) ;  • 


(5.11) 


Since  the  normal  fluid  consists  of  the  phonon  and  roton  parts,  from  Eq.  (5.11) 
we  can  obtain  directly 


.  iliOiH,!  .  _a_,  f-isl-lJ  .  [JIbL) 

4,fA(0)^^  4A(0)^  l)<A(0)j  ^  ^4^^^  L)<A(0)J 


+  .  .  .  , 


(5.13) 
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(5.14) 


We  note  that  Eqs.  (5.1) -(5. 14)  are  reduced  to  the  thermodynamic  functions  in 
Ref.  8  when  we  take  a  to  be  zero. 


VI.  First,  second  and  third  so\ind 

For  temperatures  below  0.6  K  where  the  roton  excitation  can  be 

neglected,  there  has  been  confusion  regarding  the  normal  and  anomalous 

2  7  28 

excitation  spectrums.  Some  researchers  ’  have  used  the  wrong  normal 

dispersion  to  obtain  the  temperature  variation  of  the  first  sound.  Maris  and 
29 

Massey  pointed  out  chat  the  dispersion  relation  should  be  anomalous  due  Co 
ultrasonic  attenuation. 

As  we  mentioned  earlier,  we  have  derived  the  anomalous  excitation 

15  16 

spectrums  and  introduced  a  new  approach  ’  which  includes  a  collision  term 
in  the  Boltzmann  equation,  and  obtained  the  first  and  second  sound 
simultaneously.  At  low  frequencies  such  that  «  1,  where  is  the  sound 
frequency  and  r  is  the  characteristic  time,  one  can  make  use  of  a 
hydrodynamical  approach  to  sound  propagation,  and  for  the  opposite  case  of 
»  1,  it  is  appropriate  to  make  use  of  the  kinetic  equations.  We  shall 
evaluate  the  first  and  second  sounds  in  the  kinetic  approach.  Since  we  have 
Created  the  kinetic  approach  in  detail,  we  will  briefly  review  only  the 
essential  features  and  give  the  results  on  the  basis  of  our  new  temperature - 
dependent  excitation  spectriims. 

We  solve  the  equation  of  motion  for  the  velocities  of  first  and  second 
sound  for  temperatures  below  0.6  K  where  roton  excitation  is  negligible  but 
phonon  excitation  dominates.  Our  approach  is  similar  to  that  developed  by 
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Andreev  and  Khalatnlkov  and  Dlsatnlk,  and  in  order  to  obtain  both  the 

first  and  second  sounds,  we  include  a  collision  term  and  treat  the  kinetic 

equation  for  the  phonon  distribution  function  n(p,r,t),  equation  of 

continuity,  and  equation  for  the  superfluid  velocity  v^ .  When  the  liquid  is 

slightly  perturbed  from  equilibrium,  we  assume  that  the  perturbation  terms 

from  the  distribution  mass  density  p  and  superfluid  velocity  (v^)  are 

proportional  to  exp[l(k*r  -  w^t)  ]  .  Here  k  is  the  complex  wave  vector  + 

iic^,  where  the  ratio  characterizes  the  attenuation  of  sound.  When  ^2^1 

«  1,  we  may  linearize  the  above  mentioned  three  perturbed  terms  for  n' ,  />' 

and  v'  .  To  simplify,  we  adopt  the  single  collision  model  and  use  our 

excitation  spectrum,  obtaining  two  linear  homogeneous  equations  in  p'/p  and 

V/c: 

s 


(6.1) 


where 
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*T  lcv(T)  ’  S(^j)  "  1  - 


[(w/kv(T))^  -  I]'* 


(6.3) 


A  nontrivial  solution  can  be  obtained  from  Eqs .  (6.1)  and  (6.9)  under  the 


condition  that 


-  1 

kC  ^  ^  o 


3*  -  ■  ?  f?)'  [(-■  •  rife  en 


fe  •  rife})  -  ..f?)-  0-  . 


(6.4) 


where  we  replace  the  group  velocity  by  its  thermal  average  to  obtain 


v(T)  -  C^(l  +  a'  +  37'pJ  +  ...)  . 


(6.5) 


(1  .  - 2_  .J,  ,  _  a_ 


4a2(0)  ^ 


(6.6) 


1  -  2^ 

A  ■  ”  “  "  ■  4,  ■'  -  .1 

iw  r  -  1  iw  r 

s  s 


(6.7) 


Since  the  expression  C  can  be  written  as 


^  C  +  2u^C^  ^ 

o  3  0^0  p  ,2  T  D 

o  dp 


(6.8) 
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where  -  (dP/dp)j  is  the  isothermal  sound  velocity,  we  find  the  unknown  in 
Eq .  (6.4)  to  be  w^/kC^  or  u^/kC^.  From  Eq.  (6.4)  we  can  obtain  the  solutions 
for  the  hydrodynamical  region  «  1)  and  collisionless  region  (Wg*"  »  !)• 


A.  Hydrodynamic  region  «  D 

In  the  low  frequency  limit  1)  .  can  be  approximated  as 


0  -  -1  +  2^^  -  -  lliw^r 


(6.9) 


Substitution  of  Eq  (6.9)  into  (6.4)  yields  a  quadratic  equation  in  (w^/kC)  . 

2 

Solving  Eq .  (6.4)  for  (cj^/kC)  together  with  Eq.  (6.8),  we  obtain  the  first 
and  second  sounds  and  their  attenuation  coefficients: 


^1  (o')  5  ) 

r*' -  -  1  +  -Q  4(2u+1)  +  2(2u'^+3u+l)a  +  (6u  +5u+l)^ 

Pq  I  J 


(6.10) 


(^)l  -  ^  (2u^3u.l)  .  1<2«.«} 


(6.11) 


c.d) 


P  2) 

(2a+^)  -  —  (u+1)^^  . 

) 


(6.12) 


KJ2  L 


1  +  (7u^+4u+l) 


Pj  ’ 


(6.13) 


where  a  and  $  are  given  by 


a  -  o'  +  37?^  -  55P^  , 


2  2 

rP„7  P. 


^  0-'  dp  ^  o 


(6.14) 
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with 


S  -  [2A^(0)]'^[A2  -  (1+aJ)Va^(0)]  . 


B.  CoLIlslonless  region  »  I) 


In  the  high-frequency  limit,  we  can  approximate  Eq.  (6.7)  as 


^  a.  .  -(1--)  +  2^^  -  ^ 

<p  2^  3''  2iw  r 

s 


(6.15) 


Substituting  Eq.  (6.15)  into  (6.4)  and  solving  the  quadratic  Eq .  (6.4)  for 
2 

(w^/kC)  ,  we  obtain  the  velocities  of  the  first  and  second  sounds  and 

corresponding  attenuation  coefficients: 


C  (T) 

-J; -  -  1  +  ^[(6u^+8u+3)  +  (8u^+15u+6)a  +  (10u^+llu-t-3)y91 


(6.16) 


(6.17) 


C.,(T) 


2  a2„ 

°  P. 


-  ^  [l  +  (a'+37P^-56P^)  -  ^|(5u^+12u+6)  +  ^  .  (6.18) 

L  '  o  dp  ^  ”o-* 


2  2 

Pt  ^  p^ 


(if),  -  Sif;  ['  -  2<.'-.3.p2-5»p‘)  .  i(uu2.6.2  -l“]  .  (6.19) 

l-'2  s  o  dp  o-‘ 


In  the  derivation  of  Eq .  (6.4),  we  have  adopted  the  single-collision  time 

model  for  the  collision  integral  term.  If  we  neglect  these  collisions  among 

excitations  and  solve  this  equation  for  (w  /k) ,  we  obtain 

s 
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[(u»l)2  +  1  -  -  2(utl)2<,] 

L  o  dp 


(6.20) 


We  note  that  except  for  the  a  term,  Eq.  (6.20)  is  identical  to  what  we  have 
derived  in  our  previous  paper. 

In  the  case  of  low  frequencies  (w^^  «  1).  can  apply  the  superfluid 
hydrodynamics  equations  together  with  the  dissipation  function  to  obtain  the 
first  and  second  sounds.  However,  this  derivation  is  not  related  directly  to 
our  new  temperature  variation  of  the  excitation  specti-um,  and  thus  we  will  not 
discuss  this  here. 

Since  the  third  sound  in  superfluid  helium  films  has  been  observed  in 
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both  thick  and  thin  films,  researchers  have  investigated  this  sound  in 
different  ways.  Recently  we  have  also  analyzed  thin  film  data  in  terms  of  the 
elementary  excitations  obtained  by  a  microscopic  approach.  Using  these 
previous  results  and  the  new  temperature-dependent  excitation  spectrum,  we 
will  reanalyze  the  third  sound  data. 

When  sound  waves  with  wavelength  longer  than  the  film  thickness 
propagate  along  a  helium  film,  the  normal  fluid  is  held  rigid  to  the 
substrate,  while  the  superfluid  shows  density  fluctuations.  For  the  third 
sound  velocity,  we  make  use  of  the  formula 


C?(T)  -  J  P  (T)x(T) 
J  ms 


(6.21) 


where  m  is  the  helium  mass,  p  (T)  is  the  superfluid  number  density,  and  the 

s 

adiabatic  elastic  constant  k  can  be  obtained  from  the  second  derivative  of  the 
excitation  spectrum  with  respect  to  the  surface  number  density  n.  Using  Eqs. 
(5.1),  (5.13)  and  (5.14),  we  obtain  the  third  sound  velocity 
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C^(T)  -  C^[l  +  X^(kgT)^  + 


,  -A/k  T  -A/k  T 

+  X^(kgT)  ^  e  “  +  X^CkgT)^/^  e  “  +  ... 


(6.22) 


where 


*  3 

V  -  V  -  7T  E 
o  o  10  o 


D.  -  (2jra^V*)^ 
L  o 
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D„  -  jna.  (-  £  -V  )  , 

2  4  2  o  o 


2  *  a  *  4  * 

B  «  »ra^V^[l  +  V^(3+2In(na^V^))]  . 


and  the  coefflcents  and  are  given  by 


(6.23) 


nB 

m* 


X,  - 


C(3) 

n2  3 
irD^n 


t®  ■  aojn  ]  ■  ■  Jj] 


D^n-^ 


„  3!a£i4i 

2  "  5^6 

irn 


(13.  10] 

U  ■  n2j  ’ 


S  [t  -- 


IQg 

.2 


6aD, 


2  2  2 
nDi 


1  -  ^  (1  +  nD  )(1  -  ^)1  , 

^  2D:  nD,  3 


(6.24) 
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X 


4 


3 

%  .m*  h 
n 


2r(3)q 


^5  - 


irn'^BD^ 


Ml  -  ^]  (g) 

2nDf-* 


As  in  Eq .  (4.12),  m* ,  A  and  q^  are  the  roton  effective  mass,  energy  gap  and 

minimum  point,  respectively.  These  parameters  are  given  as  a  function  of  the 

13 

potential  parameters  in  our  previous  work.  In  Eq.  (6.22)  the  first  term  is 
associated  with  the  ground-state  energy,  and  the  next  three  and  remaining 
exponential  potential  terms  are  related  to  the  phonon  and  roton  energies, 
respectively. 


VII.  Results  and  discussion 

In  the  previous  sections  we  introduced  the  eigenvalues  of  the  effective 

boson  propagator  and  the  effective  eigenvalue  for  higher  temperatures,  and 

then  we  evaluated  the  temperature -dependent  structure  factor  and  elementary 

excitation  spectrum.  We  used  these  quantities  to  obtain  thermodynamic 

functions,  fluid  density  and  various  sounds  in  thin  liquid  ^He  films. 

The  elementary  excitation  spectrum  can  generally  be  determined  from 
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neutron  scattering  measurements  and  various  sound  velocity  data.  However, 

the  shape  of  the  spectrum  in  the  very  low-momentum  region  can  not  be 

determined  well  by  such  experimental  methods,  but  it  can  be  done  indirectly 

from  specific  heat  measurements.  The  existence  of  the  anomalous  phonon 

spectrum  (concave  up)  was  revealed  in  the  first  measurements  of  the  specific 
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heat  above  0.3  K  by  Phillips  et  al.  More  recently,  Greywall  has  measured 
the  specific  heat  of  bulk  liquid  helium  with  high  precision  in  the  temperature 
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range  fron  65  to  80  mK  and  for  various  oolar  volumes  under  high  pressure  up  co 
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25  bars.  The  analysis  of  these  data  through  our  excitation  spectrum  at 
absolute  zero  temperature  gives  the  anomalous  phonon  dispersion  and  also 
agrees  with  the  results  obtained  by  sound  propagation. 

Concerning  the  specific  heat  of  thin  liquid  ^He  monolayer  films,  there 

are  several  measurements  of  the  heat  capacity  of  ^Ke  films  adsorbed  on  varioiis 
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substrates.  Among  these,  we  now  analyze  the  specific  data  of  He  monolayers 

39 

adsorbed  on  "Grafoil  graphite"  substrates  by  iising  our  Eqs.  (5.3)  and  (5.4) 
At  7(Z)  ••  0.  Here,  we  assume  that:  (1)  the  interaction  between  ^He  atoaw  is 
a  soft  potential,  (2)  the  system  is  perfectly  two  dimensional,  and  (3) 
substrate  effects  are  negligible. 

40 

Concerning  substrate  effects,  Krotscheck  considered  Chat  the  liquid 

helium  is  translationally  invariant  in  the  xy^plane,  using  three  z-dependanc 

substrate  potentials  U(z),  i.e.,  the  Aziz  potential, artificially  •weakened 
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Aziz  potential  and  Dupont-Roc  potential,  which  model  the  adhesion  force 

4 

between  the  He  atom  and  substrete.  He  then  displayed  the  one>body  densities 

piz)  for  the  above  three  potential  models  for  helium  films  of  the  different 

surface  coverages  and  made  use  of  them  to  evaluate  the  correlation  energy  and 

collective  excitation  energy.  Adopting  Krotscheck' s  method,  the  number 

density,  which  we  should  treat,  will  be  changed.  Then,  in  order  co  fit  Che 

specific  heat  data,  another  sets  of  parameters  in  Table  I  should  be  chosen. 

However,  in  the  analysis  of  the  superfluid  properties  of  the  helium  films,  we 

43 

proved  that  L/Dq  ,  where  L  is  the  healing  length  and  is  one  statistical 
atomic  layer  (D^  »  3.6A),  by  one  standard  layer  thickness  is  almost  constant 
below  1  K.  Therefore,  to  simplify  the  problem  we  have  neglected  the  substrate 


effect. 
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When  the  phonon  part  is  truncated  after  the  third  term,  we  obtain  the 
numerical  values  for  the  parameters  that  give  the  best  results  for  fitting  the 
specific  heat  data  at  various  densities,  which  are  listed  in  Table  I.  When 
higher-order  terms  of  the  phonon  part  are  included,  another  set  of  parameters 
should  be  found. 

Figure  1  illustrates  our  theoretical  specific  heat  in  comparison  with 

experimental  data.  At  temperatures  between  0.5  and  0.6  K,  the  phonon  and 

roton  contributions  are  comparable  for  the  given  densities,  and  the  latter 

becomes  dominant  as  the  temperature  Increases.  In  the  range  of  temperatures 

around  1.2  -  1.6  K,  the  phonon  part  increases  but  the  roton  part  decreases 

more  rapidly  and  hence  dominates  the  overall  behavior.  Thus  the  specific  heat 

reaches  a  maximum,  which  depends  on  the  density,  and  then  falls  sharply.  When 

we  determine  the  temperature  T^^  at  which  the  density  becomes  equal  to  the 

actual  density  of  liquid  ^He,  we  find  T^  -  1.5  K  at  the  density  of  2.79  x  10*^ 

A  Adjusting  this  value  in  accordance  with  the  ratio  of  the  theoretical  and 

experimental  values  for  bulk  ^He,  we  obtain  T^  -  1.2  K.  Bishop  and  Reppy^^ 

measured  the  superfluid  transition  temperature  of  a  thin  helium  film  adsorbed 

on  an  oscillating  substrate  and  reported  T^  -  1.215  K,  which  is  in  excellent 

2 

agreement  with  our  estimate.  In  Eq.  (5.3)  the  specific  heat  varies  as  T  at 
very  low  temperatures  when  7(Z)  is  zero.  This  variation  is  an  essential 
characteristic  of  a  two-dimensional  continuum  and  two-dimensional  Debye  model. 

Figures  2  and  3  represent  the  elementary  excitation  spectrums  and 

structure  factors,  respectively,  deduced  from  the  n\anerical  values  of  the 

parameters  at  various  densities  in  Table  I.  Because  of  the  scale  factor,  the 

upward  dispersion  is  not  shown  clearly  in  Fig.  2,  but  we  have  previously 

reported  that  the  upward  bending  becomes  stronger  as  the  density  increases  in 
17,43 

two  dimensions.  ’  In  Fig.  2,  as  the  density  increases  the  sound  velocity 
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and  roton  energy  gap  increase,  but  for  the  case  of  structure  factors  we  can 
not  find  a  general  tendency  as  for  the  excitation  spectrums. 

Let  us  now  return  to  the  temperature  variation  of  the  structure  factors. 
In  Eq.  (3.7)  the  first  bracket  expresses  the  first-order  modification  due  to 
the  finite  temperature.  The  effect  of  this  finite  temperature  comes  from  the 
effective  interaction  U^(q,^,z).  The  second  bracket  represents  the  AA^Cq) 
correction  terms.  For  low  temperatures  and  small  q  at  fixed  density,  the 
interaction  effect  of  U_j_(z,^,2)  appears  since  has  finite  growing  values 

with  increasing  temperature.  Therefore,  Eq.  (3.7)  gradually  has  larger  fixed 
values  as  q  -♦  0.  We  can  easily  confirm  that  in  the  limit  T  0  and  a(0)  ■*  0, 
Eq.  (3.7)  is  reduced  to  the  zero- temperature  expression  of  Eq.  (3.1).  We  note 
that  a(0,z)  depends  on  statistics  because  of  the  factor  7(z). 

To  draw  the  structure  factor  and  excitation  spectrum,  we  must  first 
evaluate  7(z)  as  a  function  of  temperature.  However,  7(z)  can  not  be 
represented  as  a  finite  closed  function  (see  Appendix  B)  ,  and  so  we  evaluate 
these  values  numerically  as  a  function  of  temperature.^^  Figure  4  illustrates 
7(z)  versus  temperature,  and  the  results  of  the  numerical  calculations  for 
7(z)  are  listed  in  Table  II. 

In  order  to  fit  the  specific  heat  data  at  the  density  of 
-  2  -  2 

2.79  X  10  A  ,  using  Eqs .  (5.3)  and  (5.4)  which  include  the  7(z)  term,  we 
have  chosen  the  potential  and  roton  parameters  as 

a  -  3.95  A  ,  V  -  8.75  K  ,  E  -  7.49  K  ,  q  -  0.814  ±  0.302  , 

0  ^0 

m*  -  2.04  m  ,  A  -  3.35  K  . 

Figure  5  illustrates  the  elementary  excitations  as  a  function  of  temperature 
using  the  above  parameters.  As  the  temperature  increases,  the  phonon  part 
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shows  strong  upward  bending.  However,  the  roton  energy  gap  decreases 
significantly  with  increasing  temperature,  while  the  roton  momentum 
maintains  almost  the  same  values  and  thus  seems  to  be  independent  of 
temperature . 

Figure  6  represents  the  theoretical  temperature  dependence  of  the 

structure  factor.  This  factor  takes  on  gradually  larger  values  with 

increasing  temperature  as  q  0.  The  peak  becomes  higher  for  higher 

temperatures  because  the  roton  energy  gap  decreases  with  increasing 

temperature.  This  result  agrees  with  the  temperature  dependence  of  the  liquid 

structure  function  obtained  in  a  variational  density  matrix  approach  for 
4  46 

liquid  He  at  nonzero  temperatures.  Extending  Landau's  theory,  Eendt,  Cowan 
47 

and  Yarnell  took  into  account  the  temperature  dependence  of  the  excitation 
energy  curve  in  the  temperature  range  1.1  -  1.8  K,  showing  that  the 

excitation  spectrum  generally  decreases  with  increasing  temperatures. 

48  49 

Recently  Suebka  and  Lu  adopted  a  modified  Brueckner-Sawada  method  together 

with  an  external  potential  and  reproduced  the  results  given  by  Bendt  et  al. 

Campbell  et  al^^  employed  a  variational  density  matrix  theory,  together 

with  the  minimum  principle  of  the  Helmholtz  free  energy,  to  derive  the 

elementary  excitation  spectrum  and  the  structure  factor.  Their  results  agreed 

very  well  with  the  experimentally -detexrmined  energy  of  Cowley  and  Wood^^  at 

low  momentum  (k  <  0.4  A  ^)  ,  and  their  excitation  spectmm  exhibits  strong 

anomalous  dispersions  with  increasing  temperature  at  very  short  wavelengths. 

We  note  that  Isihara  and  Samulski^  showed  that  within  the  chain  diagram 

approximation  the  excitation  spectrum  decreases,  while  the  structure  factor 

increases  with  temperature  in  bulk  liquid  He. 

We  note  that  for  the  thermodynamic  functions  the  phonon  part  contains 

both  odd  and  even  functions  of  temperature,  and  for  the  specific  heat  the 
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leading  term  is  quadratic  in  temperature,  which  is  characteristic  of  two 
dimensions.  The  roton  part  is  characterized  by  the  energy  gap  A  and  roton 
momentum  q^.  Because  of  the  exponential  factor,  it  is  small  for  low 
temperatures . 

In  Sec.  VI  we  adopted  the  single-collision  time  method  model  for  the 
collision  term  and  obtained  the  first  and  second  sounds.  Making  use  of  Eq. 
(6.8),  we  can  express  Eqs .  (6.10)  and  (6.16)  as  follows; 
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Since  the  normal  fluid  density  is  given  by  Eqs.  (5.13)  and  (5. 14)  in  both 
frequency  regions,  the  leading  term  for  the  first  sound  increases  as  T  ,  in 
striking  contrast  to  the  T*^InT  ^  increase  of  bulk  liquid  ^He .  The  absence  of 
the  logarithmic  term  in  two  dimensions  stems  from  the  angle  integrals.  Figure  7 
illustrates  the  temperature  variation  of  the  first  sound  velocity  in  both 
frequency  regions  at  very  low  and  moderately  low  temperatures.  In  the  region 
above  about  T  ~  0.9  K,  the  sound  velocity  decreases  due  to  a(^,z).  We  can 
confirm  that  the  phonon  excitation  spectrum  is  strongly  anomalous,  and  the  first 
sound  velocity  in  the  hydrodynamic  region  increases  more  than  in  the  collision 
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3  2C  2 

region.  In  this  plot,  the  parameters  u  and  (p  /C  )(d  o/dp  }  are  taken  to  be 

o  o  o 

28  29 

1.8  and  0.19,  respectively,  as  used  by  Slngh-Prakash  and  Marls.  The 
corresponding  attenuation  coefficients  in  Eqs.  (6.10)  and  (6.16)  are  given  by 

(7.3) 

(7.4) 
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The  second  sound  and  corresponding  attenuation  coefficients  given  in  Eqs. 
(6.12),  (6.13),  (6.18)  and  (6.19)  can  be  expressed  as  follows: 
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The  attenuation  coefficients  aj^(T)  and  02(T)  in  both  regions  show  the  variation 
with  P^(T)  under  the  assumption  of  constant  r.  This  result  is  similar  to  the 
However,  a^(T)  in  the  hydrodynamic  region  and  a2(T)  in  the 


bulk  case  . 
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colllslonless  region  at  temperatures  above  -0.6  K  will  depend  not  only  on 
but  also  moderately  on  a(^,Z)  and  [Eq.  (6.14)]. 

Figures  8  and  9  represent  our  theoretical  results  for  the  temperature 
variations  of  the  first  and  second  sounds  at  low  frequencies,  respectively.  We 
find  that  at  absolute  zero,  the  second  sounds  velocity  C2  is  about  1/j2  of  the 
first  sound  velocity.  As  temperature  increases  above  -0.8  K,  the  contribution 
of  a(0,Z)  to  the  first  sound  velocity  is  significant,  causing  a  dramatic 
decrease  in  the  velocity. 

We  note  in  Figure  9  that  as  temperature  increases,  the  second  sound  passes 
through  a  gentle  maximum,  reaches  a  minimum  at  about  0.6  K,  and  then  arrives  at 
another  moderate  plateau.  After  this  it  decreases  rapidly.  However,  the  effect 
of  a{0,Z)  on  the  second  sound,  which  is  different  from  that  of  the  first  sound, 
is  a  slight  increase  of  the  velocity  on  the  plateau. 

In  order  to  fit  the  data  of  the  third  sound  corresponding  to  an  atomic 
coverage  of  D  -  1.77,  we  have  adopted  the  following  potential  and  roCon 
parameters : 

a  -  2.55  A  ,  £  -  0.12  K  ,  V  -  2.17  K  ,  q  -  0.73  A  , 

000 

A  -  2 . 80  K  ,  m*  -  2 . 04  m 

These  parameters  are  slightly  different  from  those  we  have  taken  in  Ref.  17, 

which  is  due  to  a(0,z).  In  Figure  10,  we  show  that  our  microscopic  approach 

reproduces  the  temperature  variation  of  third  sound  velocity.  Our  expression 

53 

Eq.  (6.22)  is  very  similar  to  the  result  given  by  Rutledge  et  al,  who  obtained 
a  term  [see  their  Eq .  (29)]  with  a  temperature  -  dependent  coefficient. 
However,  the  coefficient  of  the  T  term  in  Eq.  (6.22)  contains  the  potential 
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parameters  as  well  as  being  temperature  dependent,  and  we  have  also  obtained 
higher -order  terms,  such  as  and  roton  terms.  Therefore,  our  derivation  is 
more  meaningful  and  should  fit  the  data  more  accurately. 

In  conclusion,  we  have  evaluated  the  temperature  variation  of  the 

4 

elementary  excitation  spectrum  of  thin  liquid  He  films  within  the  ring  diagram 
approximation.  The  effect  of  temperature  on  the  phonon  spectrxim  is  very  small 
for  very  low  temperatures,  and  thus  three-phonon  processes  do  not  play  much  of  a 
role.  However,  as  temperature  increases  from  0.6  K  to  the  vicinity  of  the  two- 
dimensional  transition  temperature  T^  -  1.21  K,  the  temperature  effect  is  very 
significant  for  the  physical  quantities  of  thin  liquid  ^He  films  as  in  bulk 
liquid  ^He. 
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Appendix  A:  Derivation  of  Eq.  (3.7) 

Substituting  Eqs.  (3.4)  and  (3.5)  into  Eq .  (3.6)  and  retaining  Just  the 
first  term,  we  have 
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where  a  -  l2y(z)/0.  We  can  rewrite  this  equation  in  two  parts  as 
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Using  the  summation  formula 
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where  -  nu[l  +  (1  +  a/nu)  ].  Through  a  similar  calculation  for  the  second 
terms  in  (A.l)  and  (A. 2)  and  making  use  of  the  summation  formula 
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we  then  obtain  Eq.  (3.7). 
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Appendix  B:  Nuaerlcal  evaltiatlon  of  't(z) 

The  total  number  of  particles  of  the  system  is  determined  by  integrating 
the  distribution  function  together  with  the  two-dimensional  density  of  states. 


ci£ 
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where  n  is  the  chemical  potential.  Expanding  the  denominator  of  Eq.  (B.l)  and 
performing  the  integration  over  energy,  we  obtain  the  simple  form 


(B.2) 
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kgT 

At  a  density  of  2.79  x  10  ^  ,  the  numerical  value  of  B  is  2.124/T,  and  we  can 
make  the  approximation  z  =•  1  -  exp[-2.12/T]  for  ^  >  1  (T  >  2.12  K) .  Therefore, 
Gj^(z)  can  be  expressed  approximately  as 

*  i 
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On  the  other  hand,  we  can  write  G^(z)  as 

^  -  j  Z  •  (B.5) 

i-1  ^ 

However,  this  function  can  not  be  expressed  as  a  finite  combination  of 

.45 

elementary  functions.  Therefore,  for  7(z)  we  must  perform  the  numerical 
calculation. 
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Table  I.  Potential  and  roton  parameters  to  fit  the  experimental  data  of 
specific  heat  in  the  microscopic  theory  presented  in  this  paper. 


nd/A^) 

a(A) 

< 

o 

m*/m 

A(k) 

0.0273 

3.565 

9.331 

8.005 

0.875  ±  0.275 

2.027 

3.29 

0.0279 

3.581 

8.369 

7.140 

0.814  ±  0.302 

2.042 

3.35 

0.0399 

3.161 

9.297 

7.943 

0.92  ±  0.350 

2.184 

4.00 

0.0419 

3.295 

8.040 

6.820 

0.90  ±  0.320 

1.891 

4.30 
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Table  II.  Numerical  values  of  7(z)  -  versus  T. 


Gi(z)  - 


201 


0.0714 

0.1732 

0.2357 

0.2857 

0.3286 

0.3895 

0.4393 

0.4507 

0.5000 

0.5571 

0.6571 

0.6893 

0.7321 

0.7750 

0.8571 

0.8940 

0.9643 

1.0914 

1.1634 

1.2120 

1.2720 

1.3790 

1.4450 

1.5150 


29.6918 
12.2401 
8 . 9944 
7.4203 
6.4516 
5.4428 
4.8258 
4.7037 
4.2400 
3.8054 
3.2262 
3.0755 
2.8957 
2.7354 
2.4734 
2.3713 
2.1984 
1.9424 
1.8222 
1.7491 
1.6666 
1.5373 
1.4671 
1.3993 


Gg(z) 

1.6449 
1 . 6448 
1.6436 
1.6398 
1.6331 
1.6169 
1.5980 
1.5930 
1.5689 
1.5369 
1.4741 
1.4528 
1.4242 
1.3955 
1.3409 
1.3169 
1.2721 
1.1955 
1.1548 
1.1286 
1.0974 
1.0451 
1.0149 
0.9844 


-  G^(z)/G^(z) 

0.0554 

0.1343 

0.1827 

0.2209 

0.2531 

0.2970 

0.3311 

0.3386 

0.3700 

0.4038 

0.4569 

0.4723 

0.4918 

0.5101 

0.5421 

0.5553 

0.5786 

0.6154 

0.6337 

0.6452 

0.6584 

0.6798 

0.6919 

0.7035 
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Figure  captions 

Figure  1.  Specific  heat  of  ^He  films  plotted  against  temperature.  The  solid 
dashed,  dotted  and  dash-dotted  lines  represent  the  theoretical 
calculations,  and  the  experimental  data  are  represented  by  +,  X,  0 
and  4. 

Figure  2.  Theoretical  excitation  spectrums  of  the  He  films  deduced  from  the 
specific  data  as  a  function  of  the  dimensionless  parameter  qa. 

Figure  3.  Structure  factor  S(q)  deduced  from  the  excitation  spectrum  as  a 
function  of  the  wave  vector  q. 

Figure  4.  Temperature  variation  of  7(2). 

Figure  5.  Temperature  variation  of  the  excitation  spectrum. 

Figure  6.  Tempertaure  variation  of  the  structure  factor  deduced  from  the 
excitation  spectrum  in  Figure  5. 

Figure  7.  First  sound  velocities  in  the  hydrodynamic  «  1) 

collisionless  regions  »  1).  The  dashed  lines  are  the 

velocities  when  a(0,z)  is  taken  into  account  in  both  regions. 

Figure  8.  First  sound  velocity  versus  temperature.  the  dashed  line  is  due 
to  the  effect  of  0(^,2). 

Figure  9.  Second  sound  velocity  versus  temperature.  The  dashed  line  is  due 
to  the  effect  of 

Figure  10.  Third  sound  velocity  for  D  -  1.77  as  a  function  of  temperature, 
where  the  circles  are  data  from  Rutledge  et  al,  and  the  solid 
curve  is  the  present  theory. 
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